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МАТЕМАТИКА 


УДК 512.622 
А.Э. МАЕВСКИЙ 


АЛГОРИТМ ПОИСКА КОРНЕЙ МНОГОЧЛЕНОВ 
С КОЭФФИЦИЕНТАМИ ИЗ КОЛЬЦА К[х у] 


Построен детерминированный алгоритм поиска корней многочленов одной пере- 
менной с коэффициентами из кольца К[х,у], где К — произвольное поле. Алгоритм 
имеет полиномиальные временную и емкостную сложности и может рассматриваться 
как распространение алгоритма Рота-Рукенштейна [2] поиска корней многочленов с ко- 
эффициентами из кольца К[х] на случай многочленов с коэффициентами из К[х, У]. 
Ключевые слова: корни многочленов, алгоритм Рота-Рукенштейна, факторизация 
многочленов, линейные делители, конечные поля. 


Введение и постановка задачи. Пусть К — поле произвольной характе- 
ристики, К[х,У/] — кольцо многочленов от переменных х, ус коэффициента- 
ми из К, К[хУ[Т| (=К[Гх, у, 7]) — кольцо многочленов от переменной Тс ко- 
эффициентами из К[Гх,У]. Под полной степенью аед(Кх,у)) многочлена Кх, у) 
(=К[х, У) будем понимать максимальную из степеней мономов, входящих в 
Кх, у), а под степенью (Т-степенью) многочлена С(х,у,Т) (ЕК х,У[7]) - мак- 
симальный показатель степени переменной 7 с которым она входит в 
О(х‚у,Т). Многочлен Кх,у) (ЕК[Гх,У]) будем называть 7-корнем многочлена 
О(х‚у,Т), если многочлен О(х,уКх,у)) нулевой. 

Рассмотрим следующую задачу: для заданного многочлена СХх,у,Т) 
(=К[х,У[Т]) и заданного целого числа & (>0) найти все 7-корни О(х,у,Т) 
полной степени не выше 4. Эта задача возникает во многих областях 
современной математики, например, в теории помехоустойчивого кодиро- 
вания при решении задач списочного декодирования [1], [2], [5]. Легко по- 
казать, что множество 

о‹а) = { Кх,) е К[х,У] | Че9(Кх,у)) < а хуКх,у)) =0} 
всех 7-корней СО(х,у, Т) полной степени не выше а находится во взаимно од- 
нозначном соответствии с множеством делителей СО(СхуТ) вида 
(Т- Кх,Уу)), где 4е9(Кх,у) < а. Поэтому исходная задача эквивалентна за- 
даче поиска всех линейных делителей многочлена О(х,у,Т) вида (Т- Кху)), 
аед(Кх,у)) < а. 

Существует несколько подходов к решению поставленной задачи. 
Например, можно использовать общие алгоритмы факторизации многочле- 
нов от нескольких переменных [3], [4], и выделить все искомые линейные 
делители специального вида. Однако вычислительная сложность при этом 
может оказаться слишком высокой, так как почти все алгоритмы фактори- 
зации многочленов от нескольких переменных вероятностные, а многочлен 
О(х,у,Т) может иметь большое количество ненужных нам линейных дели- 
телей вида (9(х,у)Т+Кх,у). В работе [5] предложен алгоритм поиска 
7-корней многочленов с коэффициентами из поля рациональных функций 
К(жи...,Хл). Так как К[ хи, ]сК(ж,....Хт) при т > 2, этот алгоритм может быть 
применен и для построения множества Фс(а). Однако он использует нетри- 
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виальную технику алгебраической геометрии и коммутативной алгебры, 
что сильно затрудняет, с одной стороны, его использование неспециалиста- 
ми, а с другой, его аппаратную или программную реализацию. 

В работе Рота и Рукенштейна [2] построен алгоритм поиска 
7-корней степени не выше & многочленов с коэффициентами из кольца 
КГ». Отметим, что в классе подобных алгоритмов алгоритм из [2] считает- 
ся одним из самых быстрых и эффективных [5]. По аналогии со схемой по- 
строения алгоритма Рота-Рукенштейна в настоящей работе построен алго- 
ритм вычисления Фо(а), обоснована его корректность и получена оценка 
его асимптотической сложности. 

Алгоритм вычисления множества ©о(а). Рассмотрим некоторый много- 
член О(х,у,Т) из кольца К[ х, У[Т]. Определим такое целое неотрицательное 
число /; что У делит О(х,у,Т), но у*! не делит О(х,у‚Т). Положим 


О(х,уТ) = Оу Т/У. 


4-1 к В . 
Пусть Кху) = а ноЛиХ у - некоторый многочлен полной 


степени & из кольца К[х,у]. Для всех целых Е [0,4] рекуррентно определим 


многочлены Их,у) (еК[х,У]), ОКхиТ) и О(хиТ) (ЕЖхУ[Т]) следующим 
образом. Положим К(х,у) = Кх,у), О(жиТ) = О (х,у,Т) = ОУ(х,у,Т), 


а 4-1 К -Е 
Иж у) = (ба(жУ) - НОУ = ох У, (1) 
О(х,иТ) = Он (х,ууТ + #1(%,0)), (2) 
К(х, и, Т) = ОЖ ТУ, (3) 


где Кл (>0)- такое целое число, что У® делит О{х,иТ), а у®*! 


Ох, у, Т). 

Лемма 1. Рассмотрим произвольное целое число Е[1,4]. Тогда 
многочлен (Т- Кх,у)) делит многочлен О’_(х, у, Т) в том и только в том слу- 
чае, когда многочлен (Т- #:(х,у)) делит многочлен Ох, у,Т). 

Доказательство. 1(>). Пусть (ТКху)) делит многочлен 
(ху, Т)=@(жу,Т)/УО. Тогда (ЕКх,У)) делит также Ох,у,Т)=@(х,ууТ+ 
+Е.(х,0)). Следовательно, 

О. (хууТ + в1(х,0)) = (Т- ху) Их) 
для некоторого (Хх,уТ) (ЕК х,У[Т]). В последнее равенство подставив 
вместо Т выражение (7- #1(х,0))/у, получим: 
О=(жуТ) = ((Т- #1(х,0)/у - ИхУ)Ижу(Т- 8х, 0))/У). 
Умножим обе части последнего равенства на у“, где / достаточно большое 
натуральное число, и, учитывая (1), получим: 
УОн(жу,Т) = (Т- #(жУ) Их УТ), ИхиТ) ЕК УТ]. 

Таким образом, (Т - #:(х,у)) делит УО„(х,у,Т), и так как у и 
(Т- Е.(ху) взаимно просты, то (Т- Е !(х,у)) делит О; (ху Т). 

2(=). Пусть (Т- в1(х,У)) делит О. (х,у, Т), то есть 

О (ху,Т) = (Т- м(хУ)ИЖуТ) 
для некоторого многочлена (Хх, у, Т) (ЕК х,У[7Т]). Используя это соотноше- 
ние, из (2) получаем: 
ОЖ Т) = (УТ + 1400) - (жу)) ИЖ УТ + 11 (%0)) = ХТ- ху) Их Т), 
где Иху/Т) = МхуУуТ + #4(х,0)). 
Следовательно, (Т - Кх,у)) делит ОжуТ), но О(жуТ) = УЗО (Ху), 
поэтому (Т- Кх,У)) делит и О’(х,у,Т). * 


не делит 
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Теорема 2. Следующие утверждения эквивалентны: 

(0)  С1- Кху) делит Чжу) 

(1) 3/в [1,9]: (Т- Ижу)) делит О" (жуТ); 

(11) — \/е [1,4]: (Т- Ижу)) делит @(жуТ); 

(У) — Тделит О-(хуТ), где О-(х у Т) = О (ху уТ+ Вы). 

Доказательство. (1)>(й). Отметим, что утверждение (Г) можно запи- 
сать как (Т- №(х,у)) делит Ох, у,Т). Тогда, согласно лемме 1, (Т- #(х,у)) 
делит О!(х,у,Т), следовательно, / в утверждении (й) можно положить рав- 
ным 1. 

(1)=(И). Следует из леммы 1. 

(1) =(м№). Пусть имеет место (1). Так как (х,у)=Ба, 
то © (ху Т) =( ТИХ, у. Т). Тогда Оз(х,у, Т)=(уТ) КхУУТ + №) и Т делит 
О (х,у,Т). 

(№)=(0. Пусть Т делит О’ж(хиТ) = О (ху уТ+ь). Используя те же 
рассуждения, что и при доказательстве первой части леммы 1, и тот факт, 
что Е(х,у) = ба, получаем, что (Т- Е(х,у)) делит Охх, у, Т). Используя далее 
лемму 1, получаем, что (Т- Кх,у)) делит О(х,у,Т). ® 

Для всех целых Е [0,4] рекуррентно определим многочлены ЛХх) (> 
К[х]), М(х, Т) (Е КХ[7Т]) следующим образом: 


Во = Их = Их“, (4) 
мик) = око.) (5) 


Отметим, что если многочлен ОХ(х,у/Т) ненулевой, то таким же будет и 
М(х,Т). 

Лемма 3. Если (7 - Кх,у)) делит О(х,у,Т), то для всех целых / [0,4] 
многочлен (7Т- РХх)) делит многочлен М(х, Г). 

Доказательство. Пусть (Т’- Кх,У)) делит О(х,у,Т). Тогда, согласно 
лемме 1, многочлен (Т - ИКх,у)) делит О%(х,‚у,Т) для всех целых Е [0,4]. 
Подставив у=0 в (Т- Кх,у)) и ОХ(х, у, Т), получим утверждение леммы. ® 

Из последней леммы вытекает следующая важная теорема. 

Теорема 4. Если (Т- Хх,у)) делит О(ху,Т), то для всех целых Е 
[0,а] коэффициенты %у...,Иа/ многочлена Кх,у) совпадают с коэффициента- 
ми №,....ВЮа- ОДНОГО Из делителей многочлена М{х,Т) вида (Т - В»), 
9е9(№(х)) = &-/.® 

Последняя теорема доставляет нам способ нахождения множества 
Ос(а). По заданному многочлену СОхуТ) вычислим многочлен 
Мхх,Т)=О(х,0,Т). Используя, например, алгоритм Рота-Рукенштейна [2], 
найдем все его 7-корни степени не выше 4. Согласно теореме 4, коэффици- 
енты Фо,...,ю любого 7-корня Кх,у) многочлена С(х,у, Т) обязательно совпа- 
дают с коэффициентами #ь,...,А. какого-либо Т-корня многочлена Мех, Т). 
Далее, для каждого полученного набора коэффициентов Аь,...,АМо согласно 
(2) и (3), вычислим многочлены О!(х,у,Т) и М(х,Т). Определив все 7-кор- 
ни ММ(х,Т) полной степени &1, получим варианты значений коэффициен- 
ТОВ №1,....я. Продолжая процесс поиска коэффициентов многочлена Хх,у), 
мы, как будет показано в теореме 5, найдем Фс(4). 

Описанный выше процесс формализуется в виде рекурсивного ал- 
горитма ЕтаКос"5. В ходе своей работы алгоритм использует две дополни- 
тельные переменные: неотрицательное целое число / определяющее уро- 
вень (глубину) рекурсии, и многочлен Хх,у) , который изменяется на каж- 
дом уровне рекурсии, а на последнем уровне становится кандидатом 
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на Т-корень. При начальном вызове алгоритма НпаКос5 оба параметра 
Ги Кх,У) должны быть нулевыми. 
Алгоритм Ет@аКВос{$ (поиск для многочлена О(х,у,Т) всех Т-корней пол- 
ной степени не выше 4): 

Вход: многочлен О(х,у,Т) (ЕК[х,У[Т]), натуральное число 4, неотри- 
цательное целое число /, многочлен Хх, у) (ЕК[х,У]); 

Выход: множество 7-корней полной степени & многочлена С(х, у, Т). 

Ш1. Найти такое целое неотрицательное число г что У делит 
О(х,у,Т), но У“! не делит С(х,у/Т); 

ШР. Положить О*(х,у,Т) := жи Т)/У; 

ШЗ. Для каждого из Ткорней Ах) степени 4- многочлена 
0(х,0, Т) выполнить: 


ШЗ.1. Положить Вр (жи Т) := О(хуТ+ВКХ)); 


ШЗ.2. Положить К(х,иТ) := В (хиут); 


ШЗ.3. Если (/= а), то 

Ш3З.3.1. Если (Т делит К(х,у,Т)), то вернуть Кх,у) + УЁКх) и выйти 
из алгоритма, иначе 

ШЗ.3.2. Выйти из алгоритма; 

Ш3.4. Выполнить ЕтаВоо($( К(х,у, Т), а, +1, Кх,у) + УВКХ). 

Конец алгоритма. 

Теорема 5. Алгоритм Рп@КВоо{, запущенный с начальными пара- 
метрами (О(х,у,Т),4,0,0), возвращает все 7-корни многочлена О(х,у,Т) сте- 
пени не выше 4. 

Доказательство. Пусть ®‹(а) (< К[х,у]) — множество всех многочле- 
нов, возвращаемых алгоритмом Нп@Кос5 (О(х,у, Т),4,0,0). Покажем, что ® 
(а) = Оса). Пусть КхУ) >е ®(а). Легко проверить, что полная степень 
Кх, у) не превышает 4. Многочлен Хх, у) построен таким образом, что Т`де- 
лит А(хиТ)=ОХ(х,ууТ+®.). Поэтому, согласно теореме 2, Кх‚у)-7-корень 
О(ху,Т) и 9 (а с Ос(а). Включение О‹(а) с ®©о(а) вытекает из теоремы 4 
и того факта, что алгоритм ЕпаКоо{ на каждом уровне рекурсии / осуще- 
ствляет полный перебор делителей вида (Т- А(х)) многочлена С”Х_х,0,Т).* 
Анализ сложности алгоритма Ет@Воо{$. На каждом уровне рекурсии / 
алгоритма ЕпаКос{5 мы находим Г-корни некоторого ненулевого многочле- 
на, и для каждого из его корней переходим на следующий уровень рекур- 
сии. Может показаться, что совокупное количество корней при переходе с 
одного уровня рекурсии на другой растет экспоненциально, однако, 
рассматриваемая ниже лемма 6 показывает, что это не так. 


Лемма 6. Пусть О(х,у,Т) = а (х, УТ (еК[х,У[Т)) - такой не- 


нулевой многочлен 7-степени БВ, что у не делит О(х,уТ), Вх) (ЕКГХ)-Т-ко- 

рень многочлена С(х,0,Т) степени не выше 4 и кратности у. Положим 

РАхуТ)=@жуут+Ы»)), Их,уТ) = РАх у, Т)/Уу, где г- такое наибольшее не- 

отрицательное целое число, что у делит Р/(х,у,Т), а У* не делит РАх,у,Т). 

Тогда Т-степень многочлена /х, Т) = Р(х,0, Т) не превосходит \. 
Доказательство. Пусть Т-степень СО(х,у, Т) равна ВБ, 


5(х,у,Т)= Оу, Т+ Их) 5, УТ", 
Р/(х,у,Т)= Оку, УТ + Щх))= "5% УУТ". 
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Поскольку Р(х) является Т-корнем кратности ‘у многочлена С(х,0,Т), то 
О(х,0,Т)=(Т - НХ)хТ) Юю длЯ некоторого Шх,Т) (еКГГТ), и 
$5$(х,0,Т)=О(х,0,Т)= (ТУШх,Т+ЕХ)). Так как Т’ делит $(х,0, Т), то 5(х,0) = 0 
при А=[0,)-1], но $,„:(х,0)=0. Последнее равнозначно тому, что у делит 
многочлены 5(х,у) при К= [0,-1] и не делит $,,.(х,0) (+0). Следовательно, 
у делит Р(х,уТ), но у! не делит Р/х‚уТ). Таким образом, ге [1,1]. 
Запишем многочлен Р(х,у, Г): 


г т: Ь р 
АжуеТ) = РАЖИТИУ = о УТ" 1+ У“ "Т". 


Теперь видно, что 7-степень Мх, Т)=Р(х,0, Т) не превышает г (= [1,/]). ® 

Следствие. Рассмотрим случай, когда на вход алгоритма ЕпаВос5 
поступает такой многочлен С(х,у,Т), что все 7-корни многочлена С(х,0, Т) 
имеют кратность 1. Тогда многочлены О%(х,0,Т), получаемые на всех по- 
следующих уровнях рекурсии, будут иметь Гстепень не выше 1, и их 
Т-корни могут быть вычислены непосредственно. ® 

Следующие две леммы являются подготовительными для теоремы 
об оценке сложности алгоритма ЕпаВос"5. 

Лемма 7. Пусть ОхуТ) (ЕК[Гх,У[Т]) — ненулевой многочлен 
степени Б. Тогда количество многочленов, возвращаемых алгоритмом 
НпаВоо(5, вызванным с параметрами (С(х,у, Т),4,0,0), не превышает ВБ, а об- 
щее количество рекурсивных обращений алгоритма Еп@Вос5 самому к 
себе не превышает Ва. 

Доказательство. Для каждого уровня рекурсии / (= [0,0]) алгоритма 
РидВос5 через ‘в, обозначим сумму Тстепеней всех многочленов 
0(х,0,Т), возникающих на шаге /1/3. Другими словами, в; равняется сумме 
всех ТГстепеней многочленов М(х,Т), возникающих в процессе поиска 
Т-корней СОхуТ). При ({=0 существует единственный многочлен 
Мох, Т)=О%(х,0,Т) Т-степени не выше В, поэтому ах < В. Согласно лемме 6, 
имеет место неравенство в, < ®;! для каждого /е [1,4]. Следовательно, ок 
Ь для каждого /е [0,4]. В итоге алгоритм Еп@Вос5 при / = а работает не 
более чем с ®о< Б многочленами, и общее количество рекурсивных вызо- 


ь 4-1 
вов РпаВоо(5 не превышает ‚0, < Ва. * 


Лемма 8. Пусть О(х,у,Т)= а (х,У)Г" (еКхУГТ) - такой 


ненулевой многочлен 7-степени В, что полная степень любого коэффици- 
ента дк(х,у) не превосходит /7. Тогда 7-степень всех многочленов на любом 
из уровней рекурсии алгоритма РНпаВоо{5 не выше В, а полная степень ко- 
эффициентов этих многочленов на уровне рекурсии / не превышает 
О(т-ЬФР). 

Доказательство. Нетрудно видеть, что шаги Ш/1, Ш/2, ШЗ.1, ШЗ.2 не 
увеличивают 7-степень. Пусть 


В(,у,Т)= Обх,у,Т+ В (х))= „а ОЬУХГ+ = УТ, 


ты ее 
В(ху,Т)= ЕЖУ, УТ)= УТ 
— многочлены, вычисляемые на шагах Ш/3.1, ШЗ.2 алгоритма ЕпаВос{5, на- 
ходящегося на Ам уровне рекурсии, и 9е9(а°” (х,у))<т(). Тогда 4е9( 


7. (ху) )<т()+Ьаед(в(х))=т()+Ба-), а Чео(у‘й,(х, у) )<т()+Ба-Н+1). Так 


11 


Раздел “Математика” 








как пК0)=т, то де9(у“й,(х, у) )\<т+ЫН-1)(а+1-]2)<т+Ба+1)(4]2+1)= = 
От + БР). 


Оценим асимптотическую сложность алгоритма Нп@Кос{5. Выбе- 
рем следующую широко распространенную в теории многочленов модель 
вычислений [2-4]. Предположим, что базовые операции поля К (сложение, 
вычитание, умножение, деление элементов), а также операции сравнения 
и присваивания имеют временную сложность С(1). Для многочленов 0(Х), 
ВХ (ЕКШ) операции сложения и вычитания имеют временную сложность 
О(тт{4е9(9(х)), аеч(#(х))}) операций поля К, а умножение имеет времен- 
ную сложность О(4ед(9(х))-деч(#х))) операций поля К. 

Теорема 9. Пусть многочлен С(х,у, Г) удовлетворяет условиям лем- 
мы 8. Тогда алгоритм НпаКоо&, вызванный с параметрами (О(х,уТ), 
а, 0, 0), имеет временную сложность О(Р@(Ь(т-Ф)? + НЬ))) операций 
поля К, где АБ) - временная сложность алгоритма поиска корней много- 
члена С(Т) (еК[Т]) степени Ь, и емкостную сложность О(ККт+Ь‹)”) эле- 
ментов поля К. 

Доказательство. На шагах Ш1, Ш2, ШЗ.2 изменения затрагивают 
только мономы, содержащие у Если коэффициенты а„(х,у) многочлена 
О(х‚у,Т) записать как элементы К[»[И|, то, согласно лемме 8, общее коли- 
чество изменяемых мономов можно оценить как О(К(тч-Ь‹Р)). Так как в 
течение работы всего алгоритма шаги ШУ, Ш2, ШЗ.2 выполняются не 6о- 
лее, чем Ба раз (лемма 7), то общий вклад этих шагов во временную слож- 
ность алгоритма равен О(Ь?@(тч+Ь‹Р)). 

Шаг Ш/3.1 удобно выполнять с помощью схемы Горнера. Нетрудно 
проверить, что в этом случае он имеет временную сложность 
ОБ«(т-Ь‹)?). Так как он выполняется в алгоритме Ра раз, его вклад в 
общую временную сложность алгоритма составляет О(6°@(т+Ьо)?). 

Для поиска всех 7-корней многочлена С\_(х,0,Т) на шаге Ш/3 можно 
воспользоваться алгоритмом Рота-Рукенштейна [2]. Его временную слож- 
ность в зависимости от В, Чи т можно оценить как О(БаВ’(т-+Ьа)+ЕЬ))), 
где АБ) — временная сложность алгоритма поиска корней многочлена СТ) 
(ЕК[7Т]) степени БВ. Следовательно, в течение работы всего алгоритма шаг 
ШЗ имеет сложность О(Р(В(т-Ьа)+НБ))). Таким образом, общая вре- 
менная сложность алгоритма ЕтКоо{5 составляет О(Б?@(Ь(т-+Ь‹)?+НЬ))) 
операций в поле К. 

В алгоритме РпаВоо(5 больше всего памяти требуется для хранения 
на каждом уровне рекурсии коэффициентов многочленов С(х,у, Т). Поэтому 
оценка О(К(т-+Ь‹)^) на емкостную сложность алгоритма следует из того, 
что общее число мономов О(х,у,Т) не превышает О((т + Ь‹)?), а общее 
число хранимых многочленов не выше В (по числу параллельно вычисляе- 
мых Т-корней). ® 

Необходимо отметить, что оценку сложности алгоритма ЕтаКос{ 
можно улучшить, с одной стороны, используя быстрые методы вычислений 
с многочленами, а с другой, более точным подсчетом числа операций. 
Выводы. В работе построен детерминированный алгоритм поиска всех 
7-корней степени не выше 4 произвольного многочлена ОхуТ) (> 
К[х,У[7]) и доказана его корректность. Показано, что алгоритм имеет по- 
линомиальные временную сложность О(Ё?а *((т-ва ?)* + АБ))) операций 
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поля К и емкостную сложность О(Кт-Ва ^)?) элементов поля К. Алгоритм 
может быть применен при решении различных задач, например, задачи 
списочного декодирования некоторых семейств алгебро-геометри-ческих 
кодов [1]. 
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А.Е. МАЕУЗКТУ 


ВООТ-ЕТМОТМС АЕСОКТТНМ РОК УМТУАВТАТЕ 
РОТУМОМТАЕ$ МТТН СОЕРРЕТСТЕМТ$ ЕВОМ К[х, у] 


Беегти с роупопта!-ите гооЯпатд аогийт Гог ипмацае ромупопта 
миИ соейаетз гогп К[х, у] мПеге К [5 а Не] ог апу сНагасей$Ис 15 сопгисе- 
е4. Оуг адогйт сап Бе меме аз ап ежепзюп оЁ Ме Ко#-Кискеп$ет'5 
гоо-Япатод а!1догийт то роупопта!5 гот К[х, УТ]. 
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